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■ The effectiveness of the Glasgow algorithm is explored via implementation in the 3d 
U(l) Gross-Neveu model and the realisation of the Goldstone mechanism in this model 
is compared and contrasted with its realisation in QCD. 

> 

^ ' 1. Introduction 

o 

\Q \ In QCD the fermion determinant is complex for chemical potential, /J non-zero, therefore 
generating an ensemble at /x is impracticable. The fermion number density, Jq, 
measures the excess of quarks relative to anti-quarks and is expected to start to rise from 
zero at some value, fi . The onset, /i D corresponds to the point where the phase of nuclear 
matter is more energetically favourable than the vacuum state (for which Jo = 0). Chiral 
symmetry restoration occurs at /x c and we expect that /Jo ~ A* c — JTip/3, where m p is 
the proton mass. Quenched simulations, on the other hand, predict /i c ~ m n /2 suggesting 
that in the limit where the bare quark mass, m q — ► chiral symmetry is restored for any 
\i 7^ 0. This result is unphysical because the pion should not couple to the baryon chemical 
potential. 

First attempts to simulate full QCD using the Glasgow method || also produced per- 
plexing results H|]. On an 8 4 lattice with bare quark mass m q = 0.01 at /3 — 5.1 we 
found fi ~ 0.1 which differs considerably from the strong coupling analysis prediction 
/i c ~ 0.65 for f3 = 5.0. Furthermore the scaling of fi D with m q was consistent with a 
Goldstone boson controlling the onset. 

This result has motivated an assessment of the effectiveness of the Glasgow method 
via its implementation in a simpler model. The method involves expansion of the Grand 
Canonical Partition Function (GCPF) as a polynomial in the fugacity variable (e M//T ). 
Consider the conventional expression for the GCPF in lattice QCD 

[dU] det (M(U, aw, m q )) e -^ (c/) (1) 

The GCPF (for fixed m q ) can be rescaled and expressed as an ensemble average of det M 
at fi = 0: 



43 



X 



= SW\ ^fco) det M ^ = °) ^ _ I det % 



/ [dU] det M(fi upd = 0) e- s °m \ det M(/x upd = 0) 



*This work was supported by EU contract no. ERBFMRXCT97-0122 



2 



where jj, up d is the chemical potential at which the statistical ensemble is updated and is 
distinct from fi meas which appears in the functional measure for the exact hybrid Monte 
Carlo (HMC) simulations and in R rw below. Note that generating the ensemble at fi up d = 
allows us to circumvent the problem of the complex action in the HMC algorithm. For 
optimum efficiency of the Glasgow method we require a large overlap between the ensemble 
generated using det M{ii upd ) and the exact ensemble generated using det M([i meas ). Let 
us define a reweighting factor R rw = d dl^(™ ea °'™) ' ^ e re ^ a ^ ve magnitude of the factor 
R rw configuration by configuration gives a measure of the overlap. If there is poor overlap 
between the simulated ensemble and the true ensemble it is conceivable that only a small 
fraction of the configurations will contribute significantly to Z (those where R rw is large 
in magnitude) in which case very high statistics would be required to extract realistic 
observables. 

The relevant features of the 3d Gross-Neveu (GN) model for \i ^ studies are that it 
has a chiral transition with a massless pion in the broken phase and it can be formulated 
such that the fermion determinant is positive definite for /i ^ 0. 

Mirroring the Glasgow reweighting technique implemented in full QCD (/i ^ 0) sim- 
ulations we performed an expansion of the Grand Canonical partition function (GCPF) 
for the 3d GN action in the fugacity variable e^ T . 

The full lattice action for the bosonized GN model with U(l) chiral symmetry is given 
in ||[7] . The functional measure used in the HMC algorithm is det(M^M). The Dirac 
fermion matrices, M and M*, can be conveniently expressed in terms of matrices G and 
V where G contains all the spacelike links while V (V') contains the forward (backward) 
timelike links 

1iM x M = Y xy + G xy + V xy e» + Vie"" ; -MA^) = + G xy + V xy e» + Vj e -" (2) 

and the term describing the Yukawa couplings of scalars to fermions is given (in terms of 
the auxiliary fields a and n on dual lattice sites x) by 

Y xy = 2i(m q + ^ ( cr ( £ ) + i£n(x)))5 xy . (3) 

° <x,x> 

The determinants of these fermion matrices are related to that of the propagator matrix 
P (following Gibbs §): 

r={- GV -v YV l) <«> 

by 

det(2zM) = e^ nt det(P - e^) ; det(2iM f ) = e^ nt det((p- 1 ) t - e - ") (5) 

Since det M has been expressed in terms of the determinant of a matrix which is diag- 
onal in e _At we can expand det M as a polynomial in e M . We can measure the averaged 
characteristic polynomial over the ensemble generated at [i = and provided that the 
coefficients are sufficiently well determined, we can use this to provide an analytic contin- 
uation to any non-zero /i. Consider a lattice with n s spatial and n t temporal dimensions. 
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Figure 1. Fermion number densities comparing simulations for three different values of 
fJ'upd with the exact hybrid Monte-Carlo simulation on a 16 3 lattice with 1/g 2 = 0.5 and 
m = 0.01. 



Determination of the eigenvalues of P nt allows us to construct the complete fugacity 
expansion for the GCPF: 
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The expansion coefficients are evaluated in the simulation and thermodynamic 

observables can be obtained from derivatives of In Z 

The simulated 3d GN U(l) model has a positive definite functional measure so we can 
choose fi up( i to investigate the influence of this choice on the observables. Exact hybrid 
Monte Carlo (HMC) simulations || showed a clear separation of the scales m n /2 and [i c 
indicating that the existence of a Goldstone mode in the spectrum of a theory need not 
precipitate chiral symmetry restoration for fi <C m p /3 in QCD. 

Does the poor overlap in the Glasgow algorithm prevent us from seeing the discontinuity 
in the number density at /i c ^> 771^/2? 

We will compare exact HMC and Glasgow method simulations on a 16 3 lattice at a four- 
fermi coupling of 1/g 2 = 0.5 and m = 0.01. We simulated the Nf = 12 (N=3 staggered) 
model. In the exact simulations for this parameter set there was a clear discontinuity 
(at fi c = 0.725(25)) in fermion number density as a function of the chemical potential. 
There was no evidence of an onset in the number density at fi ~ = 0.18(1). In fact 

ji ^ fi c for the exact HMC simulation. 

It is clear that the chemical potential fi up d at which the statistical ensemble is generated 
has a strong influence on the thermodynamic observables of the simulation. Fig. [IJ 
shows the number density for the exact HMC simulation (discussed above) and for three 
simulations with the Glasgow method: one for = 0.0 another for fi up d = 0.7( ^ fi c ) 
and finally with [i up< i = 0.8(> /i c ) The discontinuity at fi c associated with the fermion 
losing dynamical mass, which is clearly evident in the exact HMC data is not consistently 
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Figure 2. Partition function zeros for two simulations: upper diagram with [i up d = 0.0 
showing two isolated zeros emerging from lobes and lower diagram with fi up d = 0.7 showing 
an arc of zeros approaching the real axis at // ~ /i c . Simulations on 16 3 lattice with 
l/ g 2 = 0.5 and m = 0.01. 



reproduced by the Glasgow algorithm. In fact the chiral transition is seen only when 
Hupd ~ /V For fi up d > fji c , Jq reflects only the chirally symmetric phase (where m/ = m q ) 
of the exact HMC data while for [i up d = 0, Jo reflects only the phase of broken chiral 
symmetry (where mj ^> m q ). This suggests that there is insufficient overlap. 

The Lee- Yang zeros || in the complex fi plane are the zeros of Eqn. || and their 
distribution should reflect fi c . We expect the zero with the smallest imaginary part to 
approach the real axis as the lattice volume is increased. A phase transition occurs 
whenever a root approaches the real axis in the infinite volume limit. The zeros for the 
Hupd = 0.0,0.7 are plotted in Fig.^|. For \i up d = notice the two zeros emerging from 
the body of the distribution. These two isolated zeros are located at a chemical potential 
H ~ \i c . There was no evidence for \i c in J for jj, up d = so it is more likely these two 
zeros are associated with /i D rather than /i c . For fi up d = 0.7 we did see a discontinuity in 
J at /i c and in this case we see an arc of zeros forming which intersects the real //-axis 
at jj, ~ jj c . 

How do we explain the fact that [i Q ~ fi c in simulations of this model but not in QCD? 
Consider the lattice Ward identity for the chiral condensate: 



(^75^(^)^75^(3;)) = J2( tr ( G -»( x >y) G +»( x >y))) ~ {{ tr ^G(x,x))(tr^G(y,y))) 
y y 

m q 

The pion susceptibility Eqn.(0) consists of a connected channel and a disconnected 
channel. In QCD the pion has a dominant connected contribution therefore we can 
identify the pion mass in QCD with a pole in the pseudoscalar propagator, G ps , defined 
by (with G = M' 1 ): 

G ps (t) = ]T G +P (x, t)Gl p (x, t) ~ e""^ (8) 
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Figure 3. Measurements at jj, — 0.0 showing the disconnected and connected contribu- 
tions to the pion susceptibility on a 16 3 lattice with 1/g 2 = 0.5 and m = 0.01. 



Gibbs [|I|] derived a relation between the eigenvalues of P and m n in QCD. If X ps is an 
eigenvalue of P associated with the mass pole in G ps and assuming the pion susceptibility 
has a dominant connected contribution it follows that m n = 2 In | X ps | . Let us consider 
how the existence of this mass pole could affect fi Q . Notice that det M can be simply 
expressed in terms of the eigenvalues of P and Z can be similarly expressed in terms of 
its zeros aij in the e M plane: 



det M = e r 



■ n< 



4nfn t 



- on 



(9) 



Since Z = (det(M)) we see that on a single configuration oti = Aj. Note however that 
the ensemble averaged a^s are not in general the same as the ensemble averages of the 
Vs. 

The number density on a single configuration Jq ~ d In det M/dfi whereas the ensemble 
average is given by Jo ~ 91n(det M) /d[i. Presumably the unphysical singularity at fi Q ~ 
m n /2 in J associated with the pole in G ps at m 7r /2 could disappear if we achieved the ctj 
appropriate to the correct statistical ensemble and we envisage that fi Q would cancel in 
the large ensemble limit. The mass pole at m n /2 should disappear if we have a confining 
theory. 

It will be shown that in the 3d GN U(l) model the Goldstone pole forms in the discon- 
nected channel therefore the state described by G ps no longer corresponds to the Goldstone 
pion. Instead we find: 



G ~ e 



-2m f (ft=0)t 



(10) 



Since |A ps | will now correspond to the dynamical fermion mass rrif rather than m 7r /2 we 
have no reason to expect an early onset associated with the Goldstone pole in the 3d GN 
17(1) model. 

Fig. |3] shows the disconnected contribution to the Ward identity for a simulation at 
zero chemical potential. For this simulation we found (ipip) /m q ~ 40, therefore we require 
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that the sum of the connected and disconnected diagrams be of similar order so that Eqn. 
[7] is satisfied. The connected contribution was stable at a value of around 0.5 therefore 
the dominant contribution must come from the disconnected diagram. We found that the 
disconnected contribution was very noisy with large downward peaks. The data suggests 
that considering the connected contribution alone will never be sufficient to satisfy Eqn. 

We repeated our measurements for a non-zero chemical potential. We chose = 0.5 
thus ensuring that we were still in the phase of broken chiral symmetry. The results were 
consistent with those at zero chemical potential as one would expect. 

The disconnected contribution to the pion susceptibility gives a very noisy signal which 
suggests that a very long run would be required to equilibrate sufficiently to satisfy the 
lattice Ward identity. 

The Glasgow algorithm is most effective at predicting \x c for // ^ fi c when the configu- 
rations of the statistical ensemble reflect the system close to criticality . 

In QCD the pseudoscalar channel pole is formed from connected diagrams correspond- 
ing to G +ll (t)G_„{b) . A baryonic pion forms from a quark and a conjugate quark and 
condenses in this channel which explains why we find fi Q = M. This induces the 

unphysical early onset of chiral symmetry restoration in the quenched theory and will per- 
sist in the Glasgow algorithm unless either large statistics or sufficient overlap is achieved 
so that the phase of the determinant eliminates the conjugate quarks. 

In 3d GN U(l) the Goldstone mechanism is realised by a pseudoscalar channel pole 
formed from disconnected diagrams and the state G+ptyGLJt) yields a bound state of 
mass 2rrif(n = 0) which is considerably heavier than the pion. Even on individual con- 
figurations in this model we expect \i Q ~ \i c because no baryonic pion condenses 
in the connected channel. 
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